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Abstract. We study a linear spin chain which was originally introduced by Shi et al. [Phys. 
Rev. A 71 (2005), 032309, 5 pages], for which the coupling strength contains a parameter a 
and depends on the parity of the chain site. Extending the model by a second parameter (3, 
it is shown that the single fermion eigenstates of the Hamiltonian can be computed in 
explicit form. The components of these eigenvectors turn out to be Hahn polynomials with 
parameters (a, /?) and (a + l,j3— 1). The construction of the eigenvectors relies on two new 
difference equations for Hahn polynomials. The explicit knowledge of the eigenstates leads 
to a closed form expression for the correlation function of the spin chain. We also discuss 
some aspects of a g-extension of this model. 
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1 Introduction 

Consider a linear chain of N + 1 interacting fermions described by the Hamiltonian 

N-l 

H = ^ <4 {a\a k+ i + a\ +1 a k ) . (1) 

fc=0 

The lattice fermions {a k ,a\\k = 0,1,..., N} obey the common anticommutation relations, 
and Jfc expresses the coupling strength between sites k and k + 1. The Hamiltonian (JTJ) describes 
a system of N + 1 fermions on a chain with nearest-neighbour interaction (hopping between 
adjacent sites of the chain) subject to a zero background magnetic field. 

Hamiltonians of the type ([I]) appear in various contexts. In particular, spin chains of this 
type are popular as channels for short distance quantum communication, and were introduced 
by Bose [HOE]- The system then originates from a linear qubit chain with nearest neigh- 
bour interaction described by a Heisenberg XY Hamiltonian, and is being mapped into ([1]) by 
a Jordan- Wigner transformation [SJE]. In such models, the communication is achieved by state 
dynamical evolution in the spin chain, which does not require any on/off switches of the inter- 
actions between the spins, nor any modulation of external fields. Many articles dealing with 
such spin chains in the context of Bose's scheme focus on perfect transmission (or perfect state 
transfer) in these chains QUI ED H] • 

By far the most elegant and simplest scheme to realize perfect state transfer (over an arbitrary 
long chain) was proposed by Christandl et al. [7J [9] . Their choice of the modulation of the 
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coupling strengths is given by: 



Jk = JN-k-i = V(k + 1)(N - k), A; = 0,1,..., AT -1. 



(2) 



The simplicity of Christandl's model follows from the following observation. Consider first the 
single- fermion states of the system: in a single- fermion basis, the Hamiltonian H takes the 
matrix form 
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Jn-i 
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(3) 



The dynamics (time evolution) of the system is completely determined by the eigenvalues €j and 
eigenvectors ipj of this interaction matrix. It is indeed a standard technique [3 [8] to describe 
the n-fermion eigenstates of H (n < N) using the single-fermion eigenstates ipj and Slater 
determinants. In Christandl's case, determined by (|2|), the eigenvalues and eigenvectors of M 
are explicitly known. In particular, the eigenvalues are given by = — N + 2j (J = 0,1, ... , N) 
and the eigenvectors are given in terms of Krawtchouk polynomials. 

In the model of Christandl there is, for arbitrary N, an analytic (closed form) expression for 
the eigenvalues and eigenvectors of M. Such spin chains are analytically solvable [131 CE3]- It 
implies, in particular, that the correlation function at time t, 

fr,s(t) = (r\exp(-itH)\s), 

where r and s are site labels belonging to {0, 1, . . . , N} and |r), |s) denote the corresponding 
single spin states at the 'receiver' and 'sender' sites r and s, can be computed explicitly [15], HI]. 

Christandl's spin chain model allows perfect state transfer, essential for using the spin chain 
as a transmission channel. Perfect state transfer at time t = T from one end of the chain to the 
other end is expressed by \fN,o(T)\ = 1. The topic of perfect state transfer in spin chains has 
received a lot of attention [12]. Fairly easy sufficient conditions have been formulated in order 
to achieve perfect state transfer (such as mirror symmetry [5][T2]). Shi et al. showed that the 
"spectrum parity matching condition" is necessary and sufficient for perfect state transfer pQ. 
According to this condition, they found a one-parameter extension of Christandl's model in the 
case there is an even number of fermion sites in the chain; in our notation this means that A" is 
odd, i.e. 

AT = 2m + l, m £ Z+. 

The couplings in Shi's model PQ are determined by (k = 0, 1, . . . , N — 1) 



Jk 



v / (fc + l)(Af -k), if k is odd, 

y/{k + 2a + 2)(N - k + 2a + 1), if k is even. 



(4) 



Herein, a is a real parameter satisfying a > — 1 (the case of Shi actually corresponds to half- 
integer a, but for our purposes a can be any real number greater than —1). Note that for a = — h, 
Shi's model reduces to Christandl's model (i.e. reduces to ([2])), at least when A" = 2m + 1 is 
odd. In Shi et al. [I], the spectrum of the single fermion states (i.e. the eigenvalues of M with 
data determined by ([3D) was found; however no closed form expressions for the eigenvectors were 
obtained. 
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In the present paper, we will show that the eigenvectors can be expressed in terms of Hahn 
polynomials. In fact, we will first work with a two-parameter extension of Christandl's model: 

T _ j yikTi)(N-k), if^isodd, 

k \ y/{k + 2a + 2){N - k + 2ft - 1), if k is even. 

Now a and ft are real parameters satisfying a > — 1 and ft > 0. The case of Shi corresponds to 
ft = a + 1, and the case of Christandl to a = — ^, ft = \. For the general case ([5]), we obtain 
in this paper an explicit form of the eigenvalues, and an explicit form of the eigenvectors. The 
components of the eigenvectors are given by means of Hahn polynomials Q n (x; o, ft, m) I16|: 
the even components oxe proportional to Qn 

(x; a, ft, m) and the odd components to Q n (x; a + 
l,ft — l,m). In order to prove our assertions about eigenvalues and eigenvectors, we need 
some (new) difference equations for Hahn polynomials. Section [2] of this paper is devoted 
to introducing the common notation for Hahn polynomials and to proving the new difference 
equations. In Section [3] we obtain the main result of this paper: the explicit construction of the 
spectrum of M for the values ([5]) and the construction of its eigenvectors in terms of the Hahn 
polynomials. Section H] returns to the model governed by the spin chain data ([5]). Since the spin 
chain is analytically solvable, we can compute the correlation function explicitly, and determine 
under which conditions perfect state transfer is possible. Finally, in Section [5] we present the 
^-generalization of the results obtained (in terms of g-Hahn polynomials). 

Although our paper is strongly inspired by the model introduced by Shi et al., it should be 
emphasized that our results are dealing mainly with mathematical aspects of this model. In [TJ, 
the emphasis was on quantum state transfer. Our main result is to show that the eigenstates 
of the Hamiltonian ([I]) in the case of Shi, Q, or in the extended case, ([5]), can be computed in 
closed form, with coefficients given as Hahn polynomial evaluations. 

It should be mentioned that some completely different spin chain models related to Hahn 
polynomials have been considered before. The second solution of [8] is actually related to 
an interaction matrix corresponding to the Jacobi matrix of dual Hahn polynomials. In |13] . 
the interaction matrix corresponding to the Jacobi matrix of Hahn polynomials was studied, 
following some ideas of |17j . In that case, the matrix of eigenvectors U is directly a matrix of 
Hahn polynomial evaluations. However, due to the complicated coefficients in the three term 
recurrence relations, the actual coefficients in the interaction matrix become quite involved, 
see e.g. |13|. Lemma 2]. In the present paper, the main innovation comes from "doubling" the 
matrix U (hence the technique works for chains with an even length only), in a way that it 
contains Hahn polynomial evaluations of two different types (one with parameters (a, ft) and 
one with (a + 1, ft — 1)), such that the interaction matrix coefficients ([5]) are very simple. 



2 Hahn polynomials and new difference equations 

The Hahn polynomial Q n (x;a, ft,m) [T6] of degree n (n = 0, 1, . . . ,m) in the variable x, 
with parameters a > — 1 and ft > — 1, or a < —m and ft < —m is defined by |15^ I16j: 



n I a \ t? f~n,n + a + ft + l,-x \ 

Q n {x;a,ft,m)= 3 F 2 [ ;1 . (6) 

\ a + 1, — m J 

Herein, the function 3F2 is the generalized hyper geometric series |18| [19]: 



k 



a,b,c \ ^ (a)k(b)k(c)k z {7 , 

In the series is terminating because of the appearance of the negative integer — n as a nume- 
rator parameter. Note that in ([7]) we use the common notation for Pochhammer symbols [181119] 
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(a)k = a(a + l) • • • (a+ k— 1) for k = 1, 2, . . . and (a)o = 1. Hahn polynomials satisfy a (discrete) 
orthogonality relation [15] : 

m 

y] w(x; a, ft, m)Qi (x; a, P, m)Q n (x; «, P, m) = h n (a, (3,m)5i n , (8) 

x=0 

where 

/a + a;\/m + /3-a;\ 
w{x;a,f3,m) = , x = 0, l,...,m, 

\ x J \ m — x J 

, ( R _ (-1)" (n + q + g + l) m+ iQ3 + l) n n! 

nla '^' mj (2n + a + /3 + l)(a + l) n (-m) ri m! ' 

Denote the orthonormal Hahn functions as follows: 



Q n {x; a, P, m) = ===== . (9) 

y/h n {a,/3,m) 

For our construction, the essential ingredient is a set of new difference equations for Hahn 
polynomials. These relations involve Hahn polynomials of the same degree in variables x or 
x + 1, and with parameters (a, P) and (a + l,/3 — 1); in this sense it could also be appropriate 
to speak of "contiguous relations" rather than difference equations. 

Proposition 1. The Hahn polynomials satisfy the following difference equations: 

(m + (3 — x)Q n (x; a, P, m) — (m — x)Q n (x + 1; a, P, m) 

= — | Qn{x;a + 1,0 - l,m), (10 

a + 1 

(x + l)Q n (x; a + 1, P — 1, m) — (a + x + 2)Q n (x + 1; a + 1, P — 1, m) 

= -(a + l)Q n (x + l;a,P,m). (11) 

Proof. Both equations follow from a simple computation using the hyper geometric series ex- 
pression. In the case of (fit)]) , the left hand side is expanded as follows: 

(m + P — x)Q n (x; a, P, m) — (m — x)Q n (x + 1; a, P, m) 

- 1 w M ;? 1) B ;it ,) " w+m ~ * >( * — 1} - (m - - )( — 

Rewriting the expression in square brackets as [&(m — A; + 1) + (p + A;)(/c — x — 1)], the above 
sum splits in two parts: 

E (-n) fc (a + /3 + n+ l) fc (-x) fc _i -A (-re) fc (a + /3 + n + l) fc (-a;) fc 
=1 (fc-l)!(a + l) fc (-m) fc _i + ^ fc!(a + l) fc (-m) fc lP J ' 1 ' 

The first part can be brought in the following form: 

E (-n) k (a + P + n + l) fc (-x) fc _i = _ y> (-n) j+ i(a! + P + n + l) j+1 (-x)j 

n 

E 



(fc-l)!(a + l) fc (-m) fc _i ^ i!(o + l) i+ i(-m)j 

(-ra) fc (a + P + n + l) fc (-s) fc (n - k)(a + g + n + fc + 1) 
fc!(a + l) fc (-m)fc (a + fc + l) 
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So (fT2l) becomes 



E 

k=0 



(-ra)fc(a + P + n + l) k (-x) k \{n - k){a + p + n + k + I) 



k\{a + l)fe(-m) fe 



(a + fc + 1) 



+ (/3 + fc) 



E 

fc=0 



(-n)fc(a + p + n + l) fc (-x) fc r(n + a + l)(n + /?) 



fc!(a + l)fc(-m) fc 



(a + fc + l) 



(n + a + l)(ra + /3) -A (-n) fc (a + /? + n + l) fc (-s) fc 



(a + 1) 



E 

fc=0 



fc!(a + 2) fc (-m)fc 



leading to the right hand side of (fTUj) . 

The second equation is even simpler to prove. The left hand side of ([lip can be written as 



where 



k=0 V ; fc=0 



(-n)fc(a + P + n + l) fc (-a;) fc _i 



(a + k + 1)' 



(13) 



fc!(a + 2) fc _i(-m) fc 
A simple addition of the two terms in (|13p yields 

- 2^ A k {-x - 1) = -(a + 1) 2^ 

fc=0 k=0 

giving the right hand side of (|lip . 



kl(a + l) k {-m)i 



The set of difference equations will turn out to be the essential ingredient for the eigenvector 
construction in the next section. 



3 Eigenvalues and eigenvectors of the interaction matrix 

Let N = 2m + 1 be an odd integer, and consider the (N + 1) x (N + 1) interaction matrix M 
of the form ([3]) with spin chain data J k determined by ([5]), i.e. 



I y/{k + l)(2m + l-k), if k is odd, 

Jh = \ , (14) 

[ y/(k + 2a + 2)(2m + 2p-k), if is even. 

We begin with the construction of a (N + 1) x (N + 1) matrix U = (C^j)o<j,j<A r whose even 
rows are given in terms of normalized Hahn polynomials ([9]) with parameters (a, /?) and whose 
odd rows are given in terms of normalized Hahn polynomials with parameters (a + 1, P — 1). In 
order to have positive weight functions for both, we require that a > — 1 and ft > 0. 

Definition 1. The (N + 1) x (N + 1) matrix U with indices running from to N = 2m + 1 is 
defined by 

(-1)* ~ 

U 2 i,m-j = U 2 i,m+j+i = —^-Qj(i;a,P,m), (15) 
U 2 i+i, m -j = -U 2 i+i, m +j+i = — -=-Qj(i;a+ l,p - l,m), (16) 



where i, j € {0, 1, ... , m}. 
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It is easy to verify that U is an orthogonal matrix. Indeed, let us compute U T U: 

2m+l m m 



i=0 



i=0 



(U T U)jk= UjjUjk = ^ U2i,jU2i,k + ^ U2i+l,jU2i+l,k- 

i=0 ; 

For j, k G {0, . . . , m}, (fTT|) gives 

m - 

X] ^Qm-j(r,a, P,m)Q m _ k (i;a, (3, 



(17) 



i=0 



i=0 

1 1 

2 ^m— j,Tn— k i c^m—j,m—k "jfc 



using the orthogonality (jHJ) of Hahn polynomials. For j, k G {m+1, . . . , 2m+l}, the computation 
is essentially the same and (|17p gives again <5,£. For j G {0, . . . , m} and k G {m + 1, . . . , 2m + 1}, 
(fT7|) gives 



1 = 



— m— 1 



m 



i=0 



1 = 



7 > Qrn~j{i;u + !,/? - l,m)Qjfc_ m _i(i;a + 1,^-1, 



m) 



i=0 

2 ^rn—j,k—m—l 2 ^m—j,k— m— 1 0j 

and for j G {m + 1, . . . , 2m + 1} and k G {0, . . . , m}, the result is the same. So it follows that 
{U T U)jk = Sjk, or U T U = I, the identity matrix. Hence U T is the inverse of U, so UU T = I 
holds as well. 

Now we have the main proposition. 

Proposition 2. Let M be the tridiagonal (2m + 2) x (2m + 2)-matrix 



M 



where the J k are given in (|14p . and let U be the matrix determined in Definition^ Then U is 
an orthogonal matrix: 

UU T = U T U = I. 



( 





■h 







\ 




Jo 





Jl 











Jl 






















\ 









^2m 


/ 



Furthermore, the columns of U are the eigenvectors of M, i.e. 
MU = UD, 

where D is a diagonal matrix containing the eigenvalues ej of M: 

D = diag(e , ei, . . . , e 2m +i), e m -k = -2\/ (a + k + !)(/? + fc), 



(18) 



e m+fe+ i = 2 v / (a + A: + l)(/3 + /c), fc = 0, 1, . . . , 



m. 



(19) 
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Proof. The orthogonality of U has already been proved, so it remains to verify (|18p and f)19|) . 
Now 



2m+l 

(MU).. = M ikUkj = Ji-iUi-ij + JiU i+1J . (20) 



k=0 



We have to consider (|20p in four distinct cases, according to i even or odd, and to j belonging 
to {0, 1, . . . , m} or to {m + 1, m + 2, . . . , 2m + 1}. Let us consider the case that i is odd and 
j E {0, 1, . . . , m}. Then, relabelling the indices appropriately, flT3J) and (fT5|) yield: 

(MU)2i+l,m-j = J2iU2i,m~j + "^2i+1^2i+2,m-j 

= (-lyv^VCa + i + l)(m + /3 - a, ft m) 

+ (-l) i+1 v^V(* + - + l;«,/3,m) 



— (m — + l;a,/8,m)]. 

Applying (JTUJ) , this becomes 



(-1) 



V2 W — 77^ — .-mi, ^ « ^ 7-— -T^ QiW a + l,P-l,m 



i\(m — i)\hj(a, (3,m) (a + 1) 



-2y/(a + j + l)(/3 + j)U 2 i+i, m -j = e m -jU 2 i+i,m-j = {UD) 



2i+l,m—j ' 



For i odd and j £ {?n + 1, m + 2, . . . , 2m + 1}, the computation is essentially the same. For i 
even (and the two cases for j), the computation is also similar, but now the second difference 
equation (fTTT) must be used. ■ 

Note that the spectrum of M is symmetric, consisting of the values ±2-^/ (a + k + l)(/3 + k) 
(k = 0, 1, . . . , m). Furthermore, when /3 = a + 1, the spectrum consists of integers ±2(a + fc + l). 
This latter case corresponds to the model of Shi et al. pp. 



4 Some aspects of the corresponding spin chain model 

Let us consider a spin chain ([I]) with data determined by (0). The dynamics of this system 
is described by the unitary time evolution operator exp(-itH). The transition amplitude of 
a single spin excitation from site s to site r of the spin chain is given by the time-dependent 
correlation function [H [13] 

fr,s(t) = (r\exp(-itH)\s). 

But the (orthonormal) eigenvectors of H in the single fermion mode are now known and given 
by tfj = Yl Ukj\k), i.e. the columns of the matrix U constructed in (fT5|) . (fT6j) . with Htpj = 

k=0 

Mipj = Cjtfj. Using the orthogonality of the states <pj, one finds [13j : 

JV 

fr,s(t) — ^ > U r jU s jQ i . 
3=0 



8 



N.I. Stoilova and J. Van der Jeugt 



Due to the expressions (fT5j). (fT6|) . implying U r ^ m -j = (— l) r U rtm+ j + i, it is appropriate to 
write the correlation function in the following form: 

m 

f (t) — (U TJ .p- ite m-j _l_TT . TT Tp-ite m +j + l\ 

Jrsy''/ — / \ u r,m— j u s,m— j c t t - y r,m+j+l t - y s,m+j+l c ^ 

i=o 

773 

= ^C/ rim _ j C/ s , m _,(e-^™^ + (-ly+'e**™-'). (21) 

j=0 

Now it is a matter of considering the different parities for r and s. In the case they are both 
even, one finds 



/2fc,a(<) = (-l) H V w ( fc ; /8, m)w(Z; a, /3, m) 

j_g tij{a, p,m) 

In the case the first index is odd and the second even, this becomes 



= -i(-l)* + V'i"(fc;a + - l,m)to(f;a,/J,m) (22) 

~o V«-i(« + 1)P - l,m)hj(a,p,m) 

The expressions for the case even/odd and odd/odd are similar, the main message being that due 
to the analytic expressions for the eigenvectors, we obtain explicit expressions for the correlation 
function. 

Let us examine, in this context, the transition from one end of the chain (s = 0) to the final 
end of the chain (r = N = 2m + 1). Expression (|22p reduces to: 



/Wi,o(t) = -i(-l) m V(/5)m+i(a + l)m+i 

™ (2j + a + /3 + l)(-m)j sin (2^ (a + j + + j)) 
^ (j + a + /3 + l) m+1 j! " + j + !)(/? + j) 

For general a and /3, this expression cannot be simplified further. Let us now consider the 
special case that 

P = a + 1. 

Then ([23]) reduces to 

m (— \ 

/ 2m+ i,o(i) = -2i(-l) m (a + l) m+1 ^ [ - ™» sin (2t(a + j + 1)) . (24) 

j=0 ^ ^Jrn+lJ- 

This can be rewritten as 

The last sum is of hypergeometric type 2-^1, and so it can be further simplified for special values 
of t and/or a. In particular, for t = T = vr/2, one has sin (it {a. + j + 1)) = — (— sin(-7ra). 
Using then Rummer's summation formula |18|. [19] 

_ /-m,2o + 2 \ (2q + 3) ? 
2^1 o , o! _i 



2a + m + 3' / (a + 2) 
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in the right hand side of ()25|) . this expression reduces to 

f N ,o (|) = / 2m+ i,o g) = i(-l) m sin(7ra). (26) 

Note, by the way, that also for t = 2T = it one can simplify (|25p . since now sin (2ir(a + j + 1)) = 
sin(27ra). Then, using Gauss's summation formula [T8l [T9] 

-m^a + 2 ^ (m + 1), 



2a + m + 3' / (2a + m + 3) m ' 
the right hand side of ([25]) reduces to 



/jv,o(tt) = /2m+i,oOO = -2isin(27ra)(-k 



(71 



(a + l) m+ i(m + 1) 



(2a + 2) 2m +i 

Note the importance of (|26p . Indeed, keeping in mind that a > —1, we have 



7T 



7T 



113 5 



/Af,o hr = W,o k =1 for a = (27) 



2 2' 2 2' 



So there is perfect state transfer in the chain for a assuming one of these values, at time t = ir/2. 
In fact, this corresponds to the values given by Shi et al. [T]. Note that for a = — |, the spin 
chain data further reduces to that of Christandl [7j. 

The case (|27|) corresponds to 2a+l = 21 with I a nonnegative integer. As a matter of fact, this 
can still be extended slightly. Let 2a + 1 = jj^fT with both I and k nonnegative integers. Then 
for t = T' = (2fc+ l)7r/2, the factor in §25$ becomes sin(2t(a+j + l)) = sin((2/c+l)avr) = 

(—iy +k+l , and so we have the result 



7T 



2 



/jv,o (2fe + l)- = / 2m+ i, (2fc + l)- 



7T 



2 



21 

1 for 2a + 1 = — , k,l G 

2k + 1 



This case appears already in the paper of Qian et al. [20], who use the "mirror mode concurrence" 
to find this extension of Shi's result. 

As far as perfect state transfer is concerned, our extension of Shi's model by an extra para- 
meter f3 does not give rise to any new cases. In fact, just for investigating perfect state transfer, 
the mathematical machinery developed here is not necessary: the verification of the spectrum 
parity matching condition, using the entries in the interaction matrix and the spectrum itself, 
is sufficient. The main advantage of our analysis is the explicit computation of the correlation 
function. In particular, the simplicity of the expressions (|24p and (|26p . describing the transfer 
from one end of the chain to the other end, is striking. Apart from the model of Christandl [7], 
where the general correlation function is given in [131 § 2], there are no other models with such 
an elegant and simple correlation function. 



5 On the qr-generalization of the previous results 

As the classical orthogonal polynomials of hypergeometric type have a generalization in terms 
of basic hypergeometric series, i.e. a ^-generalization, one may wonder whether the present 
construction of the tridiagonal interaction matrix M can also be generalized. This is indeed the 
case: we can present a matrix M q , whose eigenvectors are given in terms of g-Hahn polynomials, 
and whose eigenvalues are symmetric and take a simple form. In order to present these results, 
let us first briefly recall some notation related to (/-series [21]. 
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For a positive real number q 1), the q-Hahn polynomial Q n (q x ;a, f3,m\q) of degree n 
(n = 0, 1, . . . , m) in q~ x is defined by I21|: 

f a~ n aBa n+1 a~ x \ 
Q n (q- X ;a,f3,m\ q )=^ 2 ( q ' m _J Q ;q,q) . (28) 

V aq,q m J 

Herein, the function 3^2 is the basic hyper geometric series \19 \ |2"T]: 

fa,b,c \ ^ (o, b, c; q) k k 
V d,e J ^ {q,d,e;q) k 

Note that in (129j) we use the common notation for g-shifted factorials and their products |21j : 

(ai, a 2 , . . . , «u; ?)jt = («i; q)k(a,2\ q)k ■ ■ ■ (<u; 

(o;g)fc = (1 - a)(l - ag) • • • (l - ag fc_1 ) and (a) = 1. 



In (|28p . the series is terminating because of the appearance of q n in the numerator. g-Hahn 
polynomials satisfy a (discrete) orthogonality relation |15j : 

m 

^w(x]a,/3,m\q)Qi(q~ x ;a,P,m\q)Q n (q~ x -,a,(3,m\q) = h n (a, /3,m\q)Si n , 

x=0 



where 

to (a; ;a,P,m \q) = -. — , _ m , (apg) , x = 0,l,...,m, 
(q,P 1 q m ;q) x 

h , o I \ _ ("/V; g) m (g, aPq m+2 ,Pq; g) n (l - a/?q)(-ag) ra (^)_ mn 
nl a, p, m|?; 9)m(a9)m(a?) a/3?> g -m. g)n(1 _ a p q 2n+l) 1 

We shall assume that < q < 1; then the weight function is positive when < a < g" 1 and 
< /3 < g _1 . Denote the orthonormal g-Hahn functions as follows: 

r, 1 -x o I \ \Jw{x\ a, /3, m|g)Q n (g _:r ; a, /3, m|g) 

Q n (g ;a,^,m|g) = - S£ === . (30) 

y'h n {q;a,fj,m) 

Just as in Section [2j the main result needed is a set of difference equations for g-Hahn 
polynomials. 

Proposition 3. q-Hahn polynomials satisfy the following q-difference equations: 

(1 - Pq m - x )Q n (q- x ; a, 0, m|g) - (l - g^Q^g"*- 1 ; a, /3, m|g) 

(l-aq" +1 )(l-/3g")g m — , , -1 n 

= : Qn(g ;«g,/3g ,mj, (31) 

1 — ag 

(1 - ? :I ' +1 )agQ n (g" :c ;ag,/3g- 1 ,m|g) - (l - aq x+2 )Q n (q~ x ~ 1 ;aq, pq' 1 ,m\q) 

= -(l-aq)Q n (q- x - 1 ;a,P,m\q). (32) 

Proof. The proof follows the same computation as in the proof of Proposition Q3 with the 
replacement of Pochhammer symbols by the corresponding (/-shifted factorials (and keeping 
track of the appropriate powers of g). ■ 
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We now come to the construction of the tridiagonal matrix M q and the matrix of eigenvec- 
tors U. The polynomials that appear here will be g-Hahn polynomials with parameters (a, f3i) 
and (aq, (5q~ l ). So in order to have positive weight functions for both sets of polynomials, we 
shall assume: 

0<g<l, 0<a<g"\ < < 1. 

As before, let N = 2m + 1, and consider the (N + 1) x (N + 1) interaction matrix M q of the 
form (|3j) with non-zero matrix elements given by: 



J 2 fc+i = 2^(1 - g fe + 1 )(l - q m ~ J 2k = 2^(1 - aq k + x )(l - pq m ~ k )q k , (33) 

where k = 0, 1, . . . , m. The (N + l) x (iV + 1) matrix U with indices running from to N = 2m + 1 
is defined similarly as in Definition [H but in terms of g-Hahn polynomials f|30j) : 

U -u 

<- J 2i,m-j — <- J 2i,m+j+l 



2i+l,m—j 



V2 

-u 2 i+l,m+j+l ~ 



-Qjiq- 



: a 



P,m\q), 



—j^-Q^q l ;aq, 



m 



(34) 
(35) 



where i, j S {0, 1, . . . , m}. By the same argument as in Section[3l the orthogonality of the matrix 
U follows from the orthogonality of the g-Hahn polynomials and the appropriate signs. 
Then the main result in the g-case reads: 

Proposition 4. Let M q be the tridiagonal (2m + 2) x (2m + 2) -matrix 

\ 



( 


Jo 





Jo 





Ji 





Ji 





V 







are 


given 


in 





J2m 
J2m J 



(1351) . Then U is 

an orthogonal matrix, UU T = U T U = I. Furthermore, the columns of U are the eigenvectors 
of M q , i.e. 

M q U = UD, 

where D is a diagonal matrix containing the eigenvalues ej of M q : 



D = diag(e ,ei, . . . ,e 2m +i), e m -k = 
2 J (1 - aq k+1 )(l - /3q k )q m - k , 



lyj (1 - aq k + x )(l - (3q k )q m - k , 



e m+k+1 = 2y(L- aq«^){l - pq")q" l - K , k = 0, 1, . . . ,m. 

The proof of this proposition is essentially the same as that of Proposition [3l and uses the 
(/-difference equations (f3T|h (f32j) . Note that, as in the ordinary case, the spectrum of M q is 
symmetric. 

Once the explicit eigenvectors of the Hamiltonian are known in the g-generalized case ()33|) . 
one can compute the correlation function, using the same expression (12 1 1) . The expressions 
become quite involved, so we give just one example here. This is in the case of transition from 
one end (s = 0) to the other end (r = 2m + 1) of the chain, and for (3 = qa: 

f 2m+ l,o(t) = i(-irq m/2 a m/2 (aq;q) m+1 



x J^(-l)''sin(2t(l 

3=0 



aq J 



(q 



m— j+1 . 



;q)j(l + aq 



(a 2 qi+ 2 ;q) m+ i(q;q) j 

The ^-generalization does not give rise to any special cases with perfect state transfer, however. 
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6 Conclusions and outlook 

We have dealt with some mathematical aspects of a spin chain model of Shi et al. [T], which 
is a one-parameter extension of the popular spin chain introduced by Christandl et al. [7]. In 
Christandl's model, the single fermion eigenvalues and eigenvectors could easily be computed, the 
eigenvectors being related to Krawtchouk polynomials. In Shi's model, with an extra parameter 
a, there was so far no known expression of the eigenvectors. In the current paper, we have 
shown that these eigenvectors can be expressed in terms of Hahn polynomials. As a matter of 
fact, we first extend Shi's model by introducing an extra parameter /3, and then construct the 
eigenvectors. In this process two types of Hahn polynomials are involved, those with parameters 
(a, f3) and those with (a + 1, f3 — 1). These Hahn polynomials are appropriately combined in 
a matrix U, yielding the eigenvectors wanted. When /3 = a + 1, the two-parameter spin chain 
reduces to that of Shi. And when a = — |, /3 = ^, the spin chain reduces to that of Christandl. 
Note, by the way, that in this last case the eigenvectors (expressed in terms of Hahn polynomials) 
indeed reduce to the known ones (expressed in terms of Krawtchouk polynomials) . This follows 
from the fact that when a = — ^ and j3 = \ the Hahn polynomials, which are 3-F2 series, reduce 
to series according to 

(s, 8 + l,-i \ (-1)* gg) (-2i-l-m- s -l \ 
3i H 3/2, -m A )--(2s + l) ( T ) 2i H ~2m- 1 ' V " 

These reductions can be obtained, e.g., from [22\ (48)]. The 2-^1 series in the right hand side 
correspond to (symmetric) Krawtchouk polynomials (with p = 1/2). 

Due to the explicit forms of the eigenvectors, the time-dependent correlation function f r s (t) 
has been computed for the spin chains under consideration. In special cases, the expression of 
the correlation function is particularly simple, see e.g. (I26j) . 

In the construction of the eigenvectors, the main relations needed are two new difference equa- 
tions for Hahn polynomials. We have also examined the (/-generalization of these results. The 
(/-extension of these difference equations is more or less straightforward. Also the construction 
of the corresponding eigenvectors in terms of (/-Hahn polynomials has been completed. 

The extension of symmetric Krawtchouk polynomials (without a parameter a) to Hahn poly- 
nomials (with parameters (a, a + 1)) may also be used in other applications. In particular, 
we hope to extend the finite oscillator models of [22], by introducing such an extra parame- 
ter. It remains to be seen, in that case, how the underlying Lie algebra is deformed, and how 
the parameter has an influence on the finite oscillator eigenstates. This topic will be treated 
elsewhere. 
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